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Abstract 


Two-dimensional  angular  optical  scattering  (TAOS)  from  clusters  of 
polystyrene  latex  spheres  is  measured  in  the  near-forward  and  near¬ 
backward  directions.  In  both  cases,  the  scattering  pattern  contains  a  rich 
and  complicated  structure  due  to  the  interaction  and  interference  of  the 
light  among  the  primary  particles.  Calculations  are  made  for  aggregates 
similar  to  those  generated  experimentally  that  demonstrate  the  rich  struc¬ 
ture  in  the  scattering  pattern.  A  comparison  of  the  experimental  and  theo¬ 
retical  TAOS  gives  good  qualitative  agreement. 
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1.  Introduction 


Recently  there  has  been  experimental  [1-3]  and  theoretical  [2,4-7]  interest 
in  the  optical  properties  of  airborne  microparticle  aggregates.  Such  aerosols 
have  been  found  to  play  important  roles  in  myriad  fields,  including  indus¬ 
trial,  environmental,  and  health  sciences  [8],  atmospheric  remote  sensing 
[9],  and  astrophysical  studies  of  scattering  by  planetary  [10]  and  interstellar 
[11]  dust  particles.  This  interest  is  primarily  piqued  by  the  many  and/or  po¬ 
tentially  harmful  materials  that  form  clusters  of  smaller  primary  particles. 
For  example,  carbon  (soot)  particles,  which  play  an  important  role  in  the 
heat  balance  of  the  earth,  are  often  found  in  fractal  aggregates  [12].  While 
the  primary  particles  may  be  in  the  Rayleigh  regime,  the  overall  cluster  is 
not,  so  it  exhibits  fundamentally  different  optical  properties.  Other  systems 
of  interest  include  dense  clusters  of  spores  that  can  pose  a  threat  as  a  bio¬ 
logical  warfare  agent  [13].  One  possible  means  for  determining  shape  and 
structure  information  about  airborne  microparticles  is  to  investigate  their 
angular  light  scattering.  Often  this  is  done  in  only  one  angular  dimension 
(6).  However,  a  growing  number  of  researchers  are  presenting  the  angular 
intensity  distribution  of  the  scattered  light  as  a  function  of  two  angles  (6,  <j>) 
[3,7,14-17].  Since  the  angular  intensity  distribution  of  nonspherical  parti¬ 
cles  (e.g.,  microparticle  aggregates,  oblate  and  prolate  spheroids,  etc)  has 
an  azimuthal  (c p )  as  well  as  polar  (0)  dependence,  observation  of  the  two- 
dimensional  angular  optical  scattering  (TAOS)  is  useful  for  determining 
shape  and  structure  information  about  such  nonspherical  particles. 

In  this  report,  we  compare  theoretical  calculations  and  experimental  results 
of  the  TAOS  from  a  three-dimensional,  close-packed  cluster  of  spherical 
primary  particles.  First,  in  section  2,  we  discuss  the  theoretical  technique  for 
calculating  the  light  scattering  from  an  aggregate.  Section  3  highlights  the 
experimental  arrangement  and  discusses  the  generation  of  the  clusters.  We 
present  the  experimental  and  theoretical  results  in  section  4.  In  section  5  we 
summarize  our  work  and  present  ideas  for  future  investigations.  Finally, 
the  appendices  outline  the  translation  addition  theorem  for  spherical  wave 
functions  and  discuss  some  relevant  aspects  for  the  computational  work. 
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2.  Theory 


In  this  section,  we  develop  a  model  that  can  be  used  to  calculate  the  differ¬ 
ential  scattering  cross  section  of  an  aggregate  composed  of  many  dielectric 
spheres.  Mie  introduced  the  exact  solution  for  the  scattering  of  a  plane  elec¬ 
tromagnetic  wave  by  a  single  isotropic  homogeneous  sphere  in  1908  [18]. 
The  vector  wave  equation  can  be  derived  by  different  approaches,  such  as 
Hertz  vectors  [19],  the  separation  of  variable  solutions  (via  the  scalar  wave 
equation)  [20],  or  T-matrix  formalism,  introduced  by  Waterman  [21].  We 
have  chosen  this  last  approach  because  it  can  be  extended  to  derive  the 
theoretical  solution  of  the  multiple  scattering  of  light  by  an  aggregate  of 
spheres.  The  model  includes  the  multiple  scattering  that  occurs  between 
the  spheres,  and  the  formulation  is  developed  to  optimize  calculations. 

2.1  The  Vector  Wave  Equation 

In  a  spherical  coordinate  system,  the  linear-independent  solution  of  the 
vector  wave  equation  (Helmholtz  equation)  in  a  source-free,  homogeneous 
medium  is  the  independent,  divergenceless  vector  spherical  wave  func¬ 
tions  \I* cmn  {kr ,  6 ,  (f. >) ; 


mn(kr,0,<f>)  =  y/jmn  <j  — ^  {kr)  P™  (cos 6)  exp  (inuf))  9 

d 


(1) 


in  {kr)  (cos  9)  exp  (i rrufi) 


^^2mn(kr,9,(l))  =  vt 


'n  (n  +  1)  . 


kr 


jn  {kr)  P™  (cos  0)  exp  (i m<j>)  \ 


36 


P™  (cos  9)  exp  (i m<j>)  6 


+ 


mi  3 
kr  sin  6  dr 


rjn  {kr)  P™  (cos  6)  exp  (i m<j>)  , 


(2) 


where  jn{kr)  is  the  set  of  spherical  Bessel  functions  of  the  first  kind  having 
complex  argument  kr;  k  =  ^  is  the  propagation  constant  for  the  medium; 
and  P™{ cos  9)  are  the  associated  Legendre  functions  of  the  first  kind  and  of 
degree  n  and  order  m,  where  n  and  m  are  integers  defined  on  the  intervals 
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1  <  n  <  oo  and  —n  <  m  <  n  .  The  time-harmonic  convention  is  exp(-i u>t). 
The  factor  7mn  =  4^+1)^+^)!  serves  to  normalize  the  spherical  wave 
functions  with  the  angular  variables  0  and  <f>.  At  large  r,  the  scattered  fields 
can  be  expressed  with  the  vectors  '&imn(kr,  0.  </>)  that  are  obtained  by  re¬ 
placement  of  the  Bessel  functions  jn{kr)  with  the  Hankel  functions  of  the 
first  kind  h£\kr)  -  jn(kr)  +  i  yn(kr).  Then,  5Re  stands  for  "the  regular  part 
of." 

Consider  a  medium  of  real  refractive  index  no,  where  a  plane  electromag¬ 
netic  wave  characterized  by  a  propagation  constant  kG  and  wavelength  A0 
impinges  on  a  dielectric  sphere  of  complex  refractive  index  m  and  propaga¬ 
tion  constant  fci .  According  to  the  Mie  theory,  the  incident  field  (Einc,  Hjnc) 
and  the  scattered  field  (Esca,  Hsca)  can  be  expanded  in  terms  of  infinite 
series  of  spherical  vector  wave  functions: 

n=oc  m=n 

Einc  =  E  E  \mn  H"  2mn]  ~  ^ ^  Q'crmn^C^  amni  (3) 

n=l  m=—n  cr,m,n 

2  n=  00  m—n 

Hinc  —  7  V  'y  ^  a lmn^6^2mn] 

and 

n— 00  m—n 

Esca  —  E  E  [/ lmn^lmn  +  /2  mn  ^2 ran]  ~~  famn^e^  <77 nm  (4) 

n= 1  m=—  n  cr,m,n 

2  n= 00  m=n 

Hsca  =  T~  £  £  [hmn'&lmn  +  /lmti^2ronl  > 

^0°  n=l  m=-n 

where  r?0  =  7^. 

2.2  T-Matrix  for  a  Single  Dielectric  Sphere 

The  T-matrix  introduced  by  Waterman  [21]  (T)  converts  the  expansion  co¬ 
efficients  of  the  incident  wave  (aamn)  impinging  on  an  arbitrarily  shaped 
particle  to  the  expansion  coefficients  of  the  radiated  scattering  wave  ( fGmn ) 
through  the  relation 

f  =  Ta.  (5) 

The  elements  of  the  T-matrix  are  independent  of  the  incident  and  scattered 
fields.  They  depend  only  on  the  intrinsic  features  of  the  scattering  particle 
(shape,  size  parameter,  and  complex  refractive  index)  and  on  its  orientation 
with  respect  to  the  coordinate  system.  For  a  spherical  particle,  the  elements 
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of  the  T-matrix  are  greatly  simplified,  and  it  is  now  well  known  that  appli¬ 
cation  of  the  boundary  condition  on  the  spherical  interface  yields 


f'lmn.lm'n' 


-/2mn,2m'n' 


_§  s  Jn  (Vo)  [kprojn  (Vo)]'  -  jn  (Vo)  [V0jw  (Vo)]' 
mm'  nn'jn  (Vo)  [Vo/>n  (Vo)]'  -  /in  ( Vo )  [VoJn  (Vo)]'  ’ 

i^Jn  (Vo)  [Vojn  (  Vo)]  -  jn  ( Vo)  [Vojn  (  Vo)] 

(  Vo)  [Vo/bi  (  Vo)]'  -  K  (Vo)  [Vojn  (fel^o)]' 

T2mn,lm'n'  =  0  > 


(6) 

(7) 

(8) 


where  the  primes  denote  derivatives  with  respect  to  the  arguments  of  the 
Bessel  functions. 


2.3  Expansion  of  the  Incident  Plane  Wave 

The  incident  plane  wave  can  be  expanded  in  terms  of  vector  spherical  wave 
functions  for  both  polarization  states,  TE  and  TM  (transverse  electric  and 
transverse  magnetic,  respectively): 


„TM 

almn 

VCmnin+J  .  a  P™ {cos  Oi)  exp  (  lm <fo)  > 

(9) 

II 

sl 

<3 

V Cmnin+3^T>™(cos  Oi)  exp  (-i mfa) , 

(10) 

nTE  - 
ulmn 

v7 Cmnin+2^fT(cos  0*)  exp  (-i m<f>i) , 

(11) 

nTE  ___ 
u2 mn  ~ 

v/Cm^in+2  .  a  .P™(cos0i)exp(  1 
sm  uj 

(12) 

where  the  factor  (mn  =  ^(^Xn+nffi *  *s  introduced  for  consistency  with 
the  normalization  chosen  for  relations  (1)  and  (2). 

2.4  Basis  of  the  Analytical  Problem 

Consider  a  cluster  of  N  isotropic,  homogeneous,  compact,  randomly  dis¬ 
persed  spheres  with  radii  a*  and  complex  refractive  indices  n\.  The  center 
of  each  sphere  Oj  is  defined  in  a  principal  coordinate  system  R0  by  a  posi¬ 
tion  vector  Ti .  The  relative  position  vector  between  two  arbitrary  spheres  i 
and  j  is  defined  by  r ^ .  Applying  the  superposition  principle  of  electromag¬ 
netic  theory,  we  find  that  the  total  external  field  Etot  is  the  summation  of 
the  incident  wave  from  the  original  source  Einc,  with  all  the  scattered  fields 
Ei  • 

-■^sca* 

N 

Etot,  —  «e¥(*o  |r  -  r0|)a°  +  £  *{k0  |r  -  r^f  V  (13) 

i=l 
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Solving  the  total  electric  field  scattered  by  an  aggregate  of  spheres  amounts 
to  evaluating  the  scattered  field  by  each  individual  sphere  in  the  cluster.  It 
is  then  necessary  to  translate  the  entire  scattered  field  on  a  common  basis 
to  get  an  exact  single  analytical  expression  of  the  total  field  scattered  by  the 
aggregate. 

2.5  Solving  the  Multiple  Scattering  Equation 

In  multiple-scattering  theory,  the  field  incident  upon  the  zth  scatterer  is  due 
to  the  contribution  of  the  field  from  the  original  source  Einc  and  the  scat¬ 
tered  fields  E|ca  from  the  other  j  particles  in  the  cluster.  This  field,  called 
the  "excitation  field"  E*xc,  has  the  formulation 

Eixc  =  E0  +  £Eia.  (14) 

j= i 

The  T-matrix  formulation  (eq  (5))  leads  to  the  relation  EgCa  =  T^1:E;,XC, 
where  is  the  isolated-scatterer  T-matrix  for  the  zth  scatterer.  Then, 
using  the  expansion  in  terms  of  spherical  vector  wave  functions  of  each 
field,  we  have 

N 

EgCa  =  T W[$te$(k0  | r  -  r0|)a°  +  £  *(k0  |r  -  r<|)/*(JV)].  (15) 

j- 1 

In  the  above  expression  (eq  (15)),  the  incident  field  is  expanded  in  the 
principal  coordinate  system  Ro(kro,  $o,  (j> o),  whereas  the  scattered  fields 
from  each  sphere  are  expressed  in  the  coordinate  systems  of  the  scatter¬ 
ed  Rj(krj,dj,<f>j).  To  solve  the  multiple  scattering  problems,  it  is  neces¬ 
sary  to  express  both  terms  in  the  zth  scatterer  coordinate  system.  Applying 
the  translational  addition  theorem  for  the  spherical  vector  wave  functions 
(eq  (A-l)  and  (A-2)  in  app  A),  and  using  the  expansion  coefficient  represen¬ 
tation  of  each  field,  we  get  from  equation  (15) 


fj(N)  =  Ti(i)  0(ijO)a0  +  ,  i  =  1  ,-,N, 


where  and  a(i,j)  are  the  translation  matrices  for  the  incident  and 

scattered  fields,  respectively.  This  last  relation  defines  a  group  of  N- 
coupled  linear  equations  whose  unknown  variables  are  the  expansion  co¬ 
efficients  of  the  scattered  field  of  each  individual  sphere.  The  solution  can 
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be  found  with  different  techniques,  such  as  the  order  of  scattering  [4],  the 
direct  inversion  matrix  [22],  or  the  iteration  method  [23].  Nevertheless,  it 
is  convenient  to  introduce  T*W,  the  N-scatterer  T-matrix  of  the  ith  scat- 
terer.  This  formulation  includes  all  the  information  about  multiple  scat¬ 
tering  effects  due  to  the  presence  of  the  N  scatterers.  It  directly  links  the 
incident  field  Eq  to  the  scattered  field  of  the  ith  scatterer  by  the  relation 

fi(A0  =  T*(JV)a<  =  T  0)a°.  (17) 

As  F(iV)  is  defined  in  the  coordinate  system  of  the  ith  scatterer,  a2  is  also. 
Then,  from  equation  (13)  and  by  application  of  the  translational  addition 
theorem  (eq  (A-3))  once  more,  the  total  scattered  field  coefficients  in  the 
principal  coordinate  system  have  the  form 

fT(Jv)  =  £(0,  i)t  *  W£(i,  0)a°.  (18) 

i=  1 


It  is  possible  to  introduce  the  global  T-matrix  of  the  cluster  Tt(n\  which 
directly  converts  the  incident  field  coefficients  to  the  total  scattered  field 
coefficients: 


fT(iV)  _  rpT(Ar)&0 


(19) 


TT(N)  =  ^/J(0  ,i)T^N)/3(i,0). 


(20) 


The  global  T-matrix  is  independent  of  the  incident  field  and  only  depends 
on  cluster  configuration;  therefore,  it  can  be  used  in  computations  for  any 
direction  and  polarization  state  of  the  incident  wave.  It  is  also  very  useful 
for  computing  orientationally  averaged  light  scattering  through  the  use  of 
the  rotational  addition  theorem. 


Substitution  of  equation  (17)  into  equation  (16)  results  in  a  new  linear  matrix 
equation,  in  which  the  unknown  variables  are  the  N- scatterer  T-matrices 
of  each  sphere.  This  system  can  be  solved  with  an  iterative  or  a  recursive 
algorithm  [5,24].  However,  for  certain  dense  aggregates,  the  iterative  algo¬ 
rithm  cannot  converge,  because  the  method  assumes  the  primary  contri¬ 
bution  is  the  incident  field.  For  this  reason,  we  use  the  recursive  T-matrix 


6 


algorithm  (RTMA)  developed  by  Chew  as  the  basis  of  our  work.  The  results 
can  be  expressed  by  these  two  equations  [24]: 


I  -  TN+1{1)  Y,  ®(n  +  1,  i)Ti{N)fi(i,  0)a(0,  N  +  1) 


j=i 


xT 


iN+i(i) 


P(N  +  l,0)  +  £  o) 

2=1 

2.*(AT+i)^(ij0)  _  fi(jV)/?(«,0)  [l  +  a(0,AT  +  l)fAr+1(JV+1)^(iV  +  1,0)]  ,  t  <  iV. 


(21) 

(22) 


The  effect  of  the  ( N  +  1)  sphere  in  the  aggregate  acts  as  a  perturbation 
for  the  whole  system.  Then,  the  individual  multiple-scattering  T-matrices 
ijii(jv-t-i)  Df  +  i)-sphere  cluster  are  evaluated  from  the  individual 

multiple-scattering  T-matrices  T*^)  of  the  Ar-sphere  cluster. 
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3.  Experiment 


The  experimental  setup,  shown  in  figure  1,  consists  of  a  particle  genera¬ 
tor,  used  to  create  the  clusters,  and  an  apparatus  to  illuminate  the  clusters 
and  measure  their  two-dimensional  angular  scattering  pattern.  Clusters  of 
polystyrene  latex  (psl)  spheres  are  generated  with  an  ink-jet  aerosol  gen¬ 
erator  (IJAG)  [25]  from  a  solution  of  psl  spheres  and  water,  loaded  into 
an  empty  ink-jet  cartridge  (Hewlett  Packard  51612A  "Blankjet"),  which  is 
mounted  on  the  IJAG's  drying  column.  The  IJAG  controller  box  sequen¬ 
tially  triggers  (either  from  1  Hz  to  2  kHz  or  singly,  on  demand)  each  of 
the  12  nozzles  in  the  cartridge  to  produce  water  droplets  (50  /am  diameter) 
containing  psl  sphere  inclusions.  The  number  of  primary  particles  within 
the  water  droplet  follows  Poisson  statistics  and  depends  on  the  concentra¬ 
tion  of  the  slurry  within  the  cartridge.  By  varying  the  concentration,  we 
can  generate  clusters  of  different  sizes.  After  being  ejected  from  the  ink¬ 
jet  cartridge,  the  droplets  are  often  accompanied  by  small  satellite  droplets. 
These  secondary  droplets  are  removed  with  a  winnow  flow,  while  the  large 
primary  droplet  is  entrained  within  an  airflow  that  carries  it  through  the 
drying  column.  The  drying  column  is  heated  to  ~105  °C,  so  that  as  the  wet 
droplets  with  inclusions  traverse  the  column,  the  water  evaporates,  leaving 
a  dry  aggregate  of  primary  particles.  At  the  end  of  the  column  is  a  tapered 
nozzle  with  an  exit  port  of  1  mm.  The  clusters  emerge  from  the  nozzle  and 
fall  into  the  scattering  volume  of  the  measuring  apparatus. 


(a)  PMT  with  632.8-nm  (p) 

interference  filter 


Nd:YAG  beam 
(532  nm) 


(c) 


Figure  1.  TAOS  measurements:  (a)  experimental  setup,  (b)  laboratory  coordinate  system,  and  (c)  illustration  of  Abbe 
sine  condition. 
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Although  the  IJAG  may  fire  periodically,  the  changing  droplet  size  and  the 
irregular  shape  of  the  aggregates  prevents  the  clusters  from  emerging  in  a 
periodic  fashion.  Due  to  the  asynchronous  passage  of  particles  through  the 
scattering  volume,  it  is  necessary  to  detect  when  a  cluster  is  in  the  correct 
location.  The  presence  of  a  falling  cluster  is  detected  by  a  continuous-wave 
He-Ne  trigger  laser  (632.8  nm)  and  a  photomultiplier  tube  (PMT)  with  a 
narrow  bandpass  filter.  As  the  aggregate  enters  the  scattering  volume,  it 
crosses  the  He-Ne  trigger  beam  and  scatters  light  toward  the  PMT.  Detec¬ 
tion  of  a  signal  with  the  PMT  initiates  the  trigger  sequence  for  the  TAOS 
measurement. 

Upon  the  laser  receiving  a  trigger  pulse,  the  second  harmonic  (A  =  532 
nm)  of  a  Q-switched  (30-ns  pulse)  Nd:YAG  illuminates  the  falling  cluster. 
The  elastic  scattered  light  is  then  collected  in  either  the  near-forward  or 
near-backward  regions  with  an  //number  =  1.2  camera  lens  used  in  the 
Abbe  sine  condition  [26].  Under  this  condition,  the  lens  is  positioned  so 
that  the  falling  particle  is  in  the  back  focal  plane  of  the  lens.  Rays  of  light 
that  scatter  off  the  particle  and  enter  the  lens  at  an  angle  7  emerge  parallel 
with  a  displacement  h,  given  by 

h  =  /  sin  7,  (23) 

where  /  is  the  focal  length  of  the  lens.  The  emerging  parallel  light  is  then 
recorded  with  an  intensified  charge-coupled  device  (ICCD)  camera  with  a 
narrow  bandpass  filter,  so  that  only  the  scattered  532-nm  light  is  recorded. 

For  our  experiments,  the  coordinate  system  of  the  laboratory  is  taken  such 
that  the  clusters  fall  in  the  —2  direction  and  the  laser  propagates  parallel  to 
the  a:-axis  (see  fig.  1(b)).  The  Abbe  sine  condition  then  relates  the  y-z  coor¬ 
dinate  of  the  ICCD  camera  to  the  0Lab-/Lab  scattering  angles  in  the  labora¬ 
tory  frame.  The  geometry  of  the  sine  condition  is  shown  in  figure  1(c).  This 
choice  of  coordinate  system  is  nonstandard  in  the  performance  of  TAOS 
calculations.  Typically,  and  as  is  done  for  our  calculations,  the  laser  propa¬ 
gates  along  the  2-axis.  In  that  case,  the  aggregates  would  be  falling  in  the 
—x  direction.  The  correspondence  between  the  calculation  and  laboratory 
frames  is  such  that  ( x ,  y,  z) caic  ■ — » (z,  -y,  x)hab-  In  order  to  make  compar¬ 
isons  between  experimental  results  and  theoretical  calculations,  it  is  useful 
to  transform  the  calculations  into  the  laboratory-frame  coordinate  system. 
This  is  done  by  comparison  of  the  two  coordinate  systems  and  use  of  the 
sine  condition  of  equation  (23).  The  resulting  transformation  is 

#Lab  =  cos-1  (sin <9Calc  cos/calc) ,  (24) 

/Lab  =  /o  +  (-1)J  sin  -1  (-  sin  <9Calc  sin  /calc) , 

where  0Lab  and  /Lab  are  the  scattering  angles  in  the  laboratory  coordinate 
frame,  0caic  and  /calc  are  the  scattering  angles  in  the  calculation  frame,  and 
/ o  is  either  0°  with  j  =  0  for  forward  scattering  or  180°  with  j  =  1  for 
backward  scattering. 
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4.  Results  and  Discussion 


To  ensure  the  proper  operation  of  our  apparatus,  TAOS  measurements  were 
performed  on  single  psl  spheres  of  two  different  sizes.  Dilute  solutions  of 
22-  and  2.29-^m-diameter  psl  spheres  were  prepared  and  loaded  into  ink¬ 
jet  cartridges.  The  TAOS  measurements  for  both  sizes  were  consistent  with 
Mie  theory,  demonstrating  the  ring  structure  associated  with  light  scatter¬ 
ing  from  a  sphere.  With  this  verification,  we  could  confidently  perform 
TAOS  measurements  on  aggregates  of  the  2.29-/jm-diameter  psl  spheres. 

Aggregates  of  2.29-/jm-diameter  psl  spheres  (refractive  index  n  =  1.59)  are 
created  with  the  IJAG.  The  slurry  is  mixed  at  0.7  mg  of  spheres  per  ml 
of  water.  In  addition,  a  small  amount  (0.01%)  of  surfactant  (Tween-20)  is 
added  to  the  water  to  help  the  mixture's  flow  properties  (necessary  for 
passage  inside  the  narrow  channels  of  the  cartridge).  The  clusters  are  char¬ 
acterized  by  an  aerodynamic  particle  sizer  (APS),  scanning  electron  mi¬ 
croscopy  (SEM),  and  optical  microscopy.  We  determine  the  average  number 
of  spheres  per  cluster  by  turning  off  the  oven  heater  and  allowing  it  to  cool. 
Wet  droplets  containing  psl  spheres  are  allowed  to  fall  through  the  column 
and  splatter  onto  a  microscope  slide.  The  water  in  each  splat  evaporates, 
leaving  a  monolayer  of  particles  that  we  count  with  an  optical  microscope 
to  determine  the  average  number  of  primary  particles  per  cluster. 

The  uniformity  of  the  droplets  can  be  inferred  from  a  count  of  the  num¬ 
ber  of  primary  particles  included  in  a  sampling  of  droplets.  In  figure  2  we 
show  the  results  of  counting  the  number  of  2.29-/im  psl  spheres  (iVsph)  in 
each  of  227  randomly  selected  splats  of  wet  droplets.  The  fraction  of  the  227 
clusters  containing  7Vsph  spheres  is  plotted  versus  Nsph .  The  most  probable 
number  of  spheres  was  found  to  be  18.  Also  plotted  is  the  Poisson  distribu¬ 
tion  probability  for  finding  Nsp^  spheres  in  a  cluster,  given  that  the  average 
Wsph  is  18.  The  two  sets  of  data  are  in  good  agreement,  as  one  would  ex¬ 
pect  if  the  original  droplets  were  monodisperse.  The  data  points  indicating 
more  than  about  40  spheres  in  a  cluster  are  probably  where  two  or  three 
droplets  hit  the  microscope  slide  at  the  same  location.  Geometric  sizing 
with  a  scanning  electron  microscope  found  this  cluster  sample  to  have  a 
mean  geometric  diameter  of  6.7  //m.  Figure  3(a)  shows  an  SEM  micrograph 
of  a  typical  psl  cluster. 

Calculations  are  performed  for  a  close-packed  cluster  containing  primary 
particles  of  the  same  size  parameter  (xpp  =  ^  s=s  13.3)  as  those  used  in  the 
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Figure  2.  Fraction  of  227 
clusters  containing  Nsp h 
spheres  plotted  versus 
iVsph.  Solid  curve  repre¬ 
sents  Poisson  distribution 
given  that  most  probable 
number  of  spheres  is  18. 


Figure  3.  Comparison  of 

(a)  SEM  of  typical  cluster 
for  which  TAOS  meas¬ 
urements  were  performed; 

(b)  Mathematica  default 
representation  of  cluster 
used  in  numerical  calcu¬ 
lations  (primary  particle 
coordinates  are  given  in 
table  1);  (c)  representation 
of  unrotated  (r  =  0°) 
theoretical  cluster  viewed 
along  z-axis  of  calculation 
frame  (i.e.,  cluster  as  seen 
from  incident  plane  wave); 
and  (d)  realization  of 
theoretical  cluster  rotated 
by  45°  and  viewed  along 
z-axis  of  calculation  frame. 


experiment.  For  the  aggregate  size  parameter  to  be  comparable  to  the  mean 
experimental  size  parameter  (xcius.  ~  40),  13  spheres  are  arranged  as  shown 
in  figure  3(b).  The  coordinates  of  each  of  the  primary  particles  are  given 
in  table  1.  Since  falling  clusters  in  the  experiment  are  randomly  oriented, 
the  calculations  are  performed  at  different  tilt  angles  (r)  so  a  reasonable 
comparison  can  be  made.  Figure  3(b)  represents  the  Mathematica  default 
view  of  the  theoretically  generated  cluster.  The  view  along  the  calculation 
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Table  1.  Coordinates  of  „  , 

spheres  used  for  cluster  — P— ere - - - — 


X 

calculation.  1 

0.000 

0.000 

0.000 

2 

1.925 

-1.189 

0.000 

3 

0.000 

1.925 

-1.189 

4 

-1.189 

0.000 

1.925 

5 

0.000 

-1.925 

-1.189 

6 

1.189 

0.000 

1.925 

7 

-1.925 

1.189 

0.000 

8 

1.189 

0.000 

-1.925 

9 

-1.925 

-1.189 

0.000 

10 

1.925 

1.189 

0.000 

11 

0.000 

-1.925 

1.189 

12 

-1.189 

0.000 

-1.925 

13 

0.000 

1.925 

1.189 

x-axis  (i.e.,  the  view  of  the  cluster  seen  from  the  illuminating  laser)  for  the 
unrotated  cluster  (r  =  0°)  is  shown  in  figure  3(c).  Rotation  of  the  cluster 
through  r  =  5°  yields  a  slightly  different  view  from  that  seen  in  figure  3(c) 
and  is  not  shown.  The  cluster  rotated  by  r  =  45°  and  viewed  from  the  illu¬ 
minating  laser  is  shown  in  figure  3(d).  The  cluster  seen  in  this  orientation 
is  remarkably  similar  to  that  shown  in  the  SEM  micrograph  of  figure  3(a). 

The  random  orientation  of  the  aggregates  emerging  from  the  IJAG  makes 
it  impossible  to  know  for  which  angles  to  calculate  the  TAOS.  Be¬ 
cause  of  this  uncertainty  we  chose  to  look  at  three  different  orientations 
(t  =  0°,  5°,  and  45°).  To  get  a  sense  of  what  the  numerical  calculations  of 
the  TAOS  look  like  in  the  laboratory  coordinate  frame,  we  used  equation 
(24)  to  transform  the  calculations.  The  resulting  transformed  forward  and 
backward  TAOS  calculations  for  the  three  different  orientations  considered 
are  shown  in  figure  4.  Due  to  the  high-intensity  forward-scattering  peak, 
the  forward-scattering  calculations  are  plotted  on  a  log  scale.  Also  shown 
are  the  regions  observed  experimentally.  In  the  forward  direction,  the  TAOS 
measurements  are  made  for  16°  <  0Caic  <  44°,  while  the  near-backward 
observations  are  made  for  120°  <  0Caic  <  148°.  These  angular  regions  cor¬ 
respond  to  the  area  bound  by  the  white  circles.  The  actual  angular  coverage 
of  the  detector  is  denoted  by  the  white  rectangle. 

As  previously  stated,  the  random  orientation  of  the  clusters  emerging  from 
the  IJAG  makes  it  impossible  to  know  the  angles  for  which  calculations 
should  be  performed.  Our  comparisons  are  limited  to  three  different  tilt 
angles  that  allow  us  to  view  the  cluster  in  three  representative  orienta¬ 
tions:  one  of  high  symmetry,  one  with  slightly  less  symmetry,  and  one 
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Figure  4.  Forward  and 
backward  TAOS  calcu¬ 
lations  transformed  into 
laboratory  coordinate 
frame  for  three  different 
cluster  orientations:  (a)  0°, 
(b)  5°,  and  (c)  45°. 


.  (b) 


(c) 


whose  symmetry  has  been  broken.  However,  the  falling  aggregates  are 
also  free  to  rotate  in  the  0Caic  direction.  Rather  than  calculate  at  a  number 
of  different  (0Caic,  4> Calc)  orientations,  we  chose  to  simply  compare  exper¬ 
iment  and  theory  by  fixing  the  tilt  angle  (t )  and  scanning  through  (f>c&\c 
within  the  bound  regions  of  figure  4.  By  sweeping  the  white  rectangle  in 
figure  4  through  the  bound  region,  we  could  choose  a  portion  of  the  calcu¬ 
lated  TAOS  pattern  to  compare  with  the  experimental  results.  Comparisons 
were  made  at  various  values  of  </>Caic  until  reasonable  correspondence  was 
found.  This  method  for  comparing  our  results  is  supported  by  calculations 
of  light  scattering  by  the  theoretical  cluster  (fig.  3)  for  parallel  and  perpen¬ 
dicular  polarizations.  These  calculations,  which  are  equivalent  to  keeping 
the  polarization  fixed  and  rotating  the  cluster  through  90°,  show  little  dif¬ 
ference  in  the  TAOS  pattern,  and  thus  give  us  confidence  that  our  compar¬ 
ison  method  is  reasonable. 
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itive  TAOS  measurements  made  in  1 
4°,  87°  <  6*Lab  <  94°).  They  demons 
s  from  near-field  interactions  and  : 
imary  particles.  Figure  6  shows  eig] 
itatively  resemble  the  experimental 
epresentations  of  different  4>c&\c  loc. 
>ns  of  figure  4  (16°  <  0caic  <  44°),  pi 


in  a  format  such  that  the  angular  coverage  within  each  plot  corresponds 
to  the  same  area  as  the  ICCD  camera.  The  best  correspondence  is  found 
for  the  cluster  tilted  by  45°.  This  is  due  to  the  lack  of  symmetry  for  this 
cluster  orientation.  For  the  unrotated  cluster,  there  exists  a  large  amount  of 
symmetry,  since  the  primary  particles  are  arranged  in  an  almost  crystalline 
fashion.  Tilting  the  cluster  by  5°  does  little  to  remove  the  symmetry  of  the 
aggregate;  however,  some  correspondence  with  the  experiment  is  found. 

For  the  near-backward  region  (120°  <  </>Lab  <  148°,  87°  <  <?Lab  <  94°), 
eight  representative  TAOS  measurements  are  shown  in  figure  7.  They, 
too,  show  the  rich  speckle  pattern;  however,  the  intensity  patches  seem 
to  exhibit  a  more  random  orientation.  The  same  type  of  comparison  be¬ 
tween  experiment  and  calculation  is  presented  for  the  backward-scattering 
case.  Figure  8  shows  the  results  of  this  comparison  with  eight  numer¬ 
ical  calculations  for  the  TAOS  in  the  bound  angular  region  in  figure  4 
(120°  <  #caic  <  148°).  A  greater  correspondence  between  experiment  and 
theory  is  found  for  the  cluster  tilted  by  45°;  however,  in  spite  of  the  high 
symmetry,  correspondence  is  also  seen  for  the  unrotated  and  5°  tilted  clus¬ 
ter  orientations. 

Although  the  comparisons  presented  in  figures  5  through  8  are  not 
quantitative,  they  do  show  some  common  features  persistent  in  both 
the  experimental  and  theoretical  results.  Most  notably,  the  comparison 
shows  a  correspondence  in  the  frequency  of  the  speckle  features,  as 
well  as  relatively  good  agreement  in  the  size  of  the  peaks.  Previously 
reported  TAOS  measurements  on  larger  (~11  ^m)  clusters  of  the  2.29-^m- 
diameter  psl  spheres  demonstrate  a  higher  frequency  speckle  pattern  with 


Figure  7.  Eight  represen¬ 
tative  TAOS  measure¬ 
ments  in  near-backward 
scattering  direction. 
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Figure  8.  Near-backward 
TAOS  calculations 
showing  qualitative 
features  found  in 
experiments. 


a  corresponding  decrease  in  the  peak  size  [3].  Quantitative  comparisons 
between  experiment  and  theory  are  not  possible  because  of  the  limited 
information  regarding  cluster  configuration  and  orientation  at  the  time 
of  the  measurement.  Despite  the  inability  to  make  a  quantitative  com¬ 
parison,  there  are  some  features  present  in  all  the  calculations  that  are 
observed  in  the  experiment.  It  is  these  features  that  should  make  it  possible 
to  determine  the  primary  particle  size,  as  well  as  overall  cluster  size.  What 
is  needed  to  better  predict  the  TAOS  is  knowledge  of  the  cluster  orientation 
and  a  more  detailed  understanding  of  the  nature  of  aggregate  formation  so 
that  we  can  more  realistically  position  the  primary  particles  of  the  cluster 
used  in  the  calculations.  The  presence  of  internal  irregularities  (e.g.,  un¬ 
evaporated  liquid  within  the  crystal)  can  also  have  a  dramatic  effect  on  the 
TAOS  [7]. 

What  we  have  presented  is  merely  one  case  demonstrating  a  qualitative 
agreement  among  the  features  found  in  the  experimentally  observed  and 
numerically  calculated  TAOS  patterns  for  an  aggregate  of  spheres.  What 
happens,  then,  if  the  size  of  the  primary  particle  is  changed?  The  effect  of 
the  primary  particles  on  the  overall  scattering  pattern  is  quite  dramatic.  The 
close  proximity  of  the  spheres  and  the  corresponding  surface  roughness  of 
the  aggregate  lead  to  near-field  effects  and  multiple  scattering  events  that 
lead  to  the  breakdown  of  the  smooth,  continuous  contours  associated  with 
scattering  from  a  single  sphere.  Although  we  can  infer,  based  on  our  knowl¬ 
edge  of  light  scattering  from  a  sphere  and  previous  experiments  on  clusters 
[3],  that  as  the  aggregate  size  changes  the  speckle  pattern  becomes  more  or 
less  dense,  the  effect  of  size  changes  in  the  primary  particles  is  less  obvious. 
One  might  expect  that  as  the  primary  particle  becomes  smaller,  the  island 
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of  intensity  would  become  broader;  and,  in  fact,  this  appears  to  be  the  case. 
Figure  9  shows  TAOS  calculations  of  the  forward-scattering  region  in  the 
laboratory  frame  for  two  unrotated  (r  =0°)  clusters  of  the  same  size,  but 
which  contain  different  size  primary  particles.  The  effect  of  primary  parti¬ 
cle  size  is  noticeable  in  the  forward-scattering  region,  where  diffraction  is 
the  dominant  component  of  the  scattering.  In  figure  9(a),  the  cluster  con¬ 
tains  2.0-/im-diameter  spheres,  while  the  cluster  in  figure  9(b)  comprises 
2.29-/xm-diameter  spheres.  Even  for  this  modest  change  in  diameter,  one 
can  clearly  see  the  effect  that  primary  particles  have  on  the  angular  inten¬ 
sity  distribution.  The  structure  of  both  clusters  is  the  same  as  that  shown 
in  figure  3(c).  This  symmetry  is  clearly  seen  in  the  diffraction  regime  of  the 
forward-scattered  light.  The  intensity  patches  in  the  TAOS  of  figure  9(a)  are 
broader,  and  there  appear  to  be  fewer  of  them.  In  contrast,  the  TAOS  in  fig¬ 
ure  9(b)  shows  sharper  peaks,  particularly  in  the  diffraction-dominated  re¬ 
gion,  as  well  as  a  higher  speckle  density  elsewhere.  The  TAOS  calculations 
depicted  in  figure  9  are  simply  meant  to  corroborate  our  results  presented 
above;  a  more  detailed  study  of  the  nature  of  light  scattering  from  clusters 
with  different  size  primary  particles  will  be  published  at  a  later  date  [27]. 


Figure  9.  Forward  TAOS 
calculations  for  two 
clusters  of  same  size 
comprising  primary 
particles  with  different 
sizes:  (a)  cluster  containing 
2.0-//m-diameter  spheres 
and  (b)  cluster  containing 
2.29-^m-diameter  spheres. 
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5.  Summary 


Although  the  results  presented  are  not  yet  quantitative,  they  do  offer  in¬ 
sight  into  the  nature  of  light  scattering  by  an  aggregate  of  spherical  primary 
particles.  The  correspondence  between  the  experimentally  observed  two- 
dimensional  angular  optical  scattering  and  the  numerically  calculated  pat¬ 
terns  is  good.  Both  experiment  and  theory  yield  the  same  qualitative  fea¬ 
tures  within  the  same  size  angular  region.  The  difficulty  in  making  quan¬ 
titative  comparisons  lies  in  unraveling  the  nature  by  which  aggregates  of 
spheres  form  and  determining  a  means  for  discerning  the  orientation  of 
the  falling  cluster  at  the  moment  of  the  TAOS  measurement.  Futhermore, 
providing  for  larger  angular  coverage  with  the  detector  would  improve  the 
ability  to  make  quantitative  comparisons.  By  far  the  easiest  of  these  lies  in 
the  large  format  detector,  capable  of  recording  the  TAOS  over  a  wide  angu¬ 
lar  range. 

The  complicated  morphology  of  airborne  microparticle  aggregates  makes 
understanding  their  light-scattering  properties  a  challenging  task.  From  an 
experimental  viewpoint,  it  is  almost  impossible  to  know  the  exact  struc¬ 
ture,  orientation,  and  number  of  primary  particles  of  the  cluster  under 
study.  Likewise,  memory  limits  and  time  constraints  often  restrict  the  size 
of  aggregates  that  can  be  studied  computationally,  and  the  cluster  struc¬ 
ture  is  somewhat  artificial  and  "crystalline."  Despite  these  difficulties,  ex¬ 
periments  and  theory  are  rapidly  finding  a  common  ground,  and  study  can 
only  improve  in  the  future. 
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Appendix  A.  Translation  Addition  Theorem 


Stein1  and  Cruzan2  first  introduced  the  translation  addition  theorem  of 
spherical  wave  functions.  This  theorem  is  formed  by  three  different  kinds 
of  relations: 

Mi mn{kri,Qi,(j>i)  =  Y  (A-l) 

fils 

+  Oij.  4>ij)$le'f!2nv(kirj,  Oj,  (f>j)  , 

=  Y  (rU .  % >  lixv{krj ,  Qj ,  4>J ) 

+  3fte  (r ij ,  Oij ,  (j>ij ) 3fte\&2 ( krj ,  6j 

If  r  <  rQl 

mn(kri,  $%->  <f>i)  =  ^  ^2j)^e^l/^(^rj’  fyify)  (A-2) 

liv 

+  B™™ {r ij  >  @ij ?  <t>ij ) 3fte ^ 2}j,v(krj ,  #j ,  </>j  )  , 

^2mn(&ri}  $ij  <£z)  =  ^j’  j ) 3fte^ lpy  ( krj ,  ,  (j>j) 

“I"  (r ij  5  ^ij  5  <t>ij  )  ^^2^1/  {krj  5  5  )  * 

If  r  >  r07 

3fteA217(^-,  (A-3) 

\IV 

$2 mn{kri,9i,<f>i)  =  Y  (rU  >  % » ^  l/w (fcr? -  6 3 >  4>j ) 

ixv 

+  3 fte  A™^  (r ^ ,  <f>ij )  ^2^  (&r j  ?  0j  5  <^j ) * 

AJ™,  3fteA™n,  £™n,  and  3fteS^n  are  translation  coefficients  that  can  be  ex¬ 
pressed  in  terms  of  Gaunt  coefficients: 

Arnn  _  (— 1)M  a  (m,  n|-p,  i/|p)  a  (n,  i/,  p)  (A-4) 

vV  p 

x  hp  ( ka )  P™~^  (cos60k)  exp[*  (m  -  p)  0ofc], 

1S.  Stein,  "Addition  theorems  for  spherical  wave  function,"  Q.  Appl.  Math.  19,  15-24 
(1961). 

zO.  R.  Cruzan,  "Translation  addition  theorems  for  spherical  vector  wave  functions,"  Q. 
Appl.  Math.  20, 33-40  (1962). 
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Bmn  _  (-1)^+1  ^o(m,»|-/i,i/|p,p-1)6(n,i/>p)  (A-5) 

VWv  p 

x  hp  ( ka )  pp-v  (cos  90k)  exp[i  (m  -  p)  cj>ok]. 

The  Gaunt  coefficients  a(m,n|— p,  u\p)  and  a  (m,  n\—p,  u\p,p  —  1)  are  ex¬ 
pressed  in  terms  of  Wigner  3 -j  symbols  (  I ,  which  are  related  to 

\0  0  0  / 

Clebsch-Gordon  coefficients: 


a(m,n\p,u\p)  =  (2p  +  1) 


I (n  +  m)\  (y  +  p,)\ (p  —  m  —  /i)! 
(n  —  m)\ (v  —  p)\(p  +  m  +  p)\ 


n  u  p  \  I  n  up 
0  0  0  /Imp  —  m  —p 


(A-6) 


a(m,n|p,^|p,g)  =  (-l)m+#' (2p  +  1) 


(n  +  m)!  (;/  +  //)!  (p  -  m  —  p)! 
(n  —  m)!  (i/  —  p)\ (p  +  m  +  /i)! 


n  v  q 
0  0  0 


n  v  p 
m  p  —m  —p 


(A-7) 


The  evaluation  of  the  translation  coefficients  with  Stein's  formalism 
leads  to  a  summation  on  a  parameter  p  by  equations  (A-4)  and 
(A-5).  It  seems  that  since  the  work  by  Peterson  and  Strom,3  there 
has  been  some  confusion  in  the  literature  about  the  correct  values 
p  takes.  Xu  recently  partially  pointed  out  the  problem.4  Our  study 
leads  to  the  following  explanation:  Respecting  the  rules  of  addition 
of  the  kinetic  moment,  the  parameter  p  should  assume  all  the  integer 
values  p  =  n  ©  u  =  (n  +  u) ,  (n  +  u  -  1) ,  (n  +  v  -  2) , ...,  \n  -  u\ .  However, 
Wigner  3 -j  symbols  vanish  if  ji+jz+jz  is  an  odd  integer.  Then  the  contribu¬ 
tion  of  the  p  summation  to  the  A™1  coefficient  value  is  null  when,  due  to  the 

(Tl  IS  T)  \ 

,  p  takes  the  values  (n  +  u  —  1) ,  (n  +  u  —  3) 

0  0  0/ 

until  |n  —  u\  +  1.  In  the  case  of  the  B""'1  coefficient,  the  Wigner  3 -j  symbol 

3B.  Peterson  and  S.  Strom,  "T  matrix  for  electromagnetic  scattering  from  an  arbitrary 
number  of  scatterers  and  representation  of  E(3),"  Phys.  Rev.  D  8,  3661-3677  (1973). 

4Y.-L.  Xu,  "Electromagnetic  scattering  by  an  aggregate  of  spheres:  Far  field,"  Appl.  Opt. 
36,  9496-9508  (1997). 


is  \  n  v  q  I,  where  q  =  p  -  1.  Then  the  contribution  of  the  p  summa- 

\0  0  o) 

tion  is  null  when  p  takes  on  the  following  values:  (n  +  v) ,  (n  +  v  -  2)  until 
\n  —  v\.  That  is,  in  the  evaluation  of  the  A™  coefficients,  p  takes  on  the  val¬ 
ues  from  | n  -  u\  to  (n  +  v)  in  steps  of  2;  and  in  the  evaluation  of  the  B™1 
coefficients,  p  takes  on  the  values  from  \n  —  v\  +  1  to  (n  +  u)  in  steps  of  2. 

The  index  n  is  linked  to  the  multipole  expansion  representative  of  each 
dielectric  sphere  of  the  cluster.  The  computation  of  the  T-matrix  through 
the  extended  boundary  condition  (EBC)  technique  requires  a  truncation  of 
an  infinite-dimensional  integral  matrix  equation  to  a  value  n  =  nmax.  This 
value  must  be  large  enough  to  correctly  represent  the  physically  relevant 
partial  waves,  but  not  too  large  to  be  useful  numerically.  It  is  usually  ad¬ 
mitted  that  Wiscombe's  criterion  allows  an  accurate  determination  of  this 
value.5  There  are  as  many  different  nmax  as  the  number  of  different  types 
and  sizes  of  spheres.  We  truncate  the  incident  and  all  the  series  expansions 
of  the  scattered  fields  at  the  larger  nmax  value.  Then,  /  and  o  are  column 
vectors  of  dimension  nmax  x  1. 


5W.  J.  Wiscombe,  "Improved  Mie  scattering  algorithms,"  Appl  Opt.  19, 1505-1509  (1980). 
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Appendix  B.  Truncation  of  the  Translation  Addition  Theorem  Series 


The  computation  of  the  N -spheres  problem  assumes  that  the  series  ex¬ 
pansions  of  the  translation  addition  theorem  (eq  (A-l),  (A-2),  and  (A-3) 
in  app  A)  are  uniformly  convergent.  Therefore,  we  can  truncate  them  at 
v  =  i/max,  assuming  that  the  resulting  error  is  small  enough  if  vm3X  is  suffi¬ 
ciently  large. 

The  a(i,  j )  matrices  describe  the  continuity  of  the  tangential  components  of 
all  the  scattered  electric  fields  across  each  surface  of  the  ith  sphere.  Then,  in 
an  iterative  formalism,  pmax  represents  the  approximation's  order  of  mul¬ 
tiple  scattering  effects  in  the  system.  Its  maximal  value  for  convergence 
is  approximately  equal  to  nmax.  In  the  recursive  algorithm,  all  multiple¬ 
scattering  effects  are  taken  into  account  during  the  matrix  inversion  of  re¬ 
lation  (21)  in  the  main  report.  Following  the  analytical  coherence  of  the 
recursive  T-matrix  algorithm  (RTMA)  formalism,  vmsx  =  nmax.  Therefore, 
a(i,j)  matrices'  dimensions  are  [2nmax(nmax  +  2)]2. 

Relation  (17)  can  also  be  described  as  a  phase  shift  of  the  incident  plane 
wave,  when  it  travels  from  origins  00  to  Op. 

a1  =  a0  exp  [ihnc  •  r*] .  (B-l) 

Starting  with  i/max  =  nmax,  we  calculate  the  absolute  variation  between  the 
numerical  results  of  equation  (B-l)  and  equation  (17)  in  the  main  report 
for  each  vector's  component,  due  to  the  variation  of  umax.  The  convergence 
speed  also  depends  on  the  distance  and  orientation  between  the  coordinate 
systems.  Nevertheless,  at  umax  =  nmax,  the  maximal  absolute  variation  is 
approximately  1.  When  vmax  =  2nmax,  it  is  only  6  to  10,  and  it  is  not  until 
vmax  =  3nmax  that  an  absolute  variation  of  10  to  15  occurs  in  the  best  cases. 
The  dimension  of  the  /? (i,  0)  matrices  is  4nmax(nmax  +  2)  x  vmsx{vmax  +  2). 

Because  of  the  analytical  coherence  of  the  RTMA  formalism,1  the  series 
(eq  (A-3))  has  to  be  truncated  at  the  same  order  umgx  as  equation  (A-2). 
The  dimension  of  the  matrices  /?(0,  *)  is  4^max(^max  +  2)  x  nmax(«max  +  2). 
The  limits  of  convergence  that  we  have  encountered  for  the  (3(i,  0)  matrices 
are  still  present  and  are  propagated  into  numerical  evaluations. 

1w.  c.  Chew,  Waves  and  Field  in  Inhomogeneous  Media,  IEEE  Press  Series  on  Electromag¬ 
netic  Waves,  1990. 
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We  can  see  that  if  the  indices  nmax  and  ^max  become  too  large,  the  numeri¬ 
cal  evaluations  of  the  matrices  a(i,  j),  j 3(i ,  0),  and  /?((),  i)  are  limited  by  the 
free  memory  available  on  the  computer,  realistic  computing  time,  and  nu¬ 
merical  convergence  of  the  series.  Taking  these  three  factors  into  account, 
we  see  that  the  computation  of  the  total  scattered  field  by  a  random  cluster 
in  terms  of  a  single  expansion  on  one  principal  coordinate  system  is  im¬ 
practicable  most  of  the  time.  Consequently,  the  scattering  parameter  of  the 
sphere  cluster  (cross  sections  and  scattering  matrix)  cannot  be  evaluated. 

To  circumvent  this  difficulty,  Borghese  et  al,2  Mackowski  and  Mishchenko,3 
and  Fuller4  derived  explicit  expressions  for  the  total  cross  sections  (scatter¬ 
ing,  extinction,  and  absorption)  of  the  sphere  cluster,  without  any  need  for  a 
single  representation  of  the  total  scattered  field.  This  formalism  is  called  the 
sphere-centered  formulation.  Moreover,  Xu5  recently  introduced  a  formal¬ 
ism  to  obtain  a  single  expression  of  the  total  scattered  field  in  the  far-field 
region.  This  approach  is  the  most  suitable  to  obtain  differential  scattering 
cross  sections  of  an  aggregate,  without  taking  the  convergence  of  the  trans¬ 
lation  addition  theorem  series  (eq  (A-3))  into  account.  Indeed,  for  an  arbi¬ 
trary  point  M  in  the  far-field  region  at  the  vector  position  r  in  the  principal 
system  and  r,  in  the  coordinate  system  of  the  ith  sphere,  one  can  consider 
that  these  two  vectors  are  parallel.  By  introducing  A  *  as  the  path  variation 
due  to  the  relative  position  between  the  principal  and  the  ith  systems,  we 
obtain  r*  =  r  —  A*  with  A,  =  er  x  ro* ,  where  er  is  a  unit  vector  pointed 
in  the  direction  of  M.  Then  it  can  be  shown  that  in  the  far-field  region,  the 
translation  coefficient  can  be  written  as 

(kroi,$oi,  cf>oi)  =  exp  (  ifcAj) ,  (B-2) 

MeB™(krQl,doi,M  =  0. 

This  formulation  leads  to  a  simplification  of  the  analytical  formulation  of 
the  /3(0,  i)  matrices.  Finally,  the  total  scattered  field  by  the  entire  cluster  of 
equation  (18)  can  be  written  in  the  far-field  region  as 

fT(iV)  _  ^  exp  (— ifcAj)  fi(Ar).  (B-3) 

i~  1 


2F.  Borghese,  P.  Denti,  R.  Saija,  G.  Toscano,  and  O.  I.  Sindoni,  "Multiple  electromagnetic 
scattering  from  a  cluster  of  spheres.  I.  Theory/'  Aerosol  Sci.  Technol.  3,  227-235  (1984). 

3D.  W.  Mackowski  and  M.  I.  Mishchenko,  "Calculation  of  the  T  matrix  and  the  scattering 
matrix  for  ensembles  of  spheres,"  J.  Opt.  Soc.  Am.  A 13, 2266-2278  (1996). 

4K.  A.  Fuller,  "Scattering  and  absorption  cross  sections  of  compounded  spheres. 
I.  Theory  for  external  aggregation,"  J.  Opt.  Soc.  Am.  A 11, 3251-3260  (1994). 

5Y.-L.  Xu,  "Electromagnetic  scattering  by  an  aggregate  of  spheres:  Far  field,"  Appl.  Opt. 
36,  9496-9508  (1997). 


Appendix  C.  Recursive  T-Matrix  Formulation  in  the  Far-Field  Region 


rj,N+l(N+l) 


rju(jV-f-l) 


Because  of  the  problems  introduced  in  appendices  A  and  B,  we  introduce 
a  new  model  of  multiple  scattering  computation  by  an  aggregate  of  spheres 
based  on  the  recursive  T-matrix  algorithm  (RTMA)  algorithm  developed  by 
Chew,1  modified  to  bypass  the  convergence  problem  of  equation  (A-l)  in 
appendix  A.  Indeed,  the  evaluation  of  the  incident  field  on  each  sphere 
is  made  through  a  phase  shift  formalism  instead  of  a  matrix  calculation. 
Then,  in  the  second  stage,  to  circumvent  the  convergence  problem  of 
equation  (A-3),  we  use  the  far-field  approximation  introduced  by  Xu.2  The 
fundamental  equations  of  this  model  are 


I  -  tN+1{1)  Y,  +  1,  i)Ti{N)a(i,  N  +  1) 


i—  1 


-1 


(C-l) 


ij.JV+1(1) 


a(N  +  1,  i)T1^  exp  [ikinc  (rjv+i  -  r^)] 


1=1 


_  f-j-ii(jv)  jj  ^  jy  +  1)T N+1(N+1 )  exp  Jik.nc  (rjV+1  _  r,)]  ,  i  <  N,  (C-2) 


N+ 1  . 

f(JV+!)  =  ^  fi(W)  exp  (-ifeAi) , 


2  =  1 

p(N+i)  _  -jn(JV-i-i)ao exp  (ikinc  x  r,) . 


(C-3) 

(C-4) 


'w.  C.  Chew,  Waves  and  Field  in  Inhomogeneous  Media,  IEEE  Press  Series  on  Electromag¬ 
netic  Waves,  1990. 

2Y.-L.  Xu,  "Electromagnetic  scattering  by  an  aggregate  of  spheres:  Far  field,"  Appl.  Opt. 
36,9496-9508(199 7). 
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